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| We study spin excitations in a two dimensional p-wave superconductor with A = d(ki ± ife) 

symmetry in the context of the newly discovered superconducting S^RuCU. The polarization and 
spectrum of spin wave excitations are identified and their experimental consequences are discussed. 
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I. INTRODUCTION 

o 

Recent discovery of superconducting S^RuCU generates a lot of efforts to determine the pairing symmetry of 
a possible unconventional superconducting order parameter in this system. Being a 4d-orbital analogue Q of the 
Qh' high T c cuprate superconductors, Sr2Ru04 has the same layered perovskite structure as La2CuC>4 and becomes 
superconducting below T c « 1.5K f^H^]- Despite the structural similarity, it behaves very differently from the copper 
oxides. The normal state below 50K can be well described by a quasi-2D Landau-Fermi liquid: The resistivity shows 
T 2 behavior in the a-b plane and c-direction with a large anisotropy ratio jlj . Quantum oscillations [[| revealed three 
cylindrical Fermi surface sheets in accordance with the band structure calculations ||. There are three 4d-orbitals 
(d X y,d yz ,d zx ) of the Ru 4+ -ions which form three bands crossing the Fermi level. There are two electron-like and one 
hole- like Fermi surfaces ||. On the other hand, the mass enhancement (« 4) in this material is not small indicating 
the existence of strong correlation The facts that it has a large mass enhancement and a related material, SrRuC>3, 
shows ferromagnetism lead to the proposal that the superconducting state of S^RuC^ is formed by odd-parity pairing 
(spin-triplet pairing) presumably in the p-wave channel |?]-|9| . 
' There have been a number of experiments indicating unconventional superconductivity in Sr2Ru04. The transition 
temperature was found to be very sensitive to non-magnetic impuritie s JlCfl and nuclear quadrupole resonance (NQR) 
' found no Hebel-Slichter peak pT[ . Even though these experiments p0|-|l2| suggested that the pairing symmetry is non- 
s-wave, they could not determine the pairing symmetry itself. More recently 17 0-Knight shift in NMR was measured 
and it is consistent with the spin triplet superconductivity jl3j with d parallel to the c axis as the case of supcrfluid 
3 He -A |Q. Here d is the unit vector of the triplet order parameter fji| ]. On the other hand, the specific heat data 
are consistent with usual s-wave superconductors if we subtract a persistent T-linear term p!^JIq |. In the earlier 
ON ■ experiment the ratio of this coefficient of the T-linear term to the one in the normal state 70/7./V was larger than 0.5. 
This initiated the proposal of possible nonunitary state in Sr2Ru04 However more recent data show that 

70/ 'jN < 0.25. This means the nonunitary state is untenable. The specific heat data further imply that the energy 
gap should be almost independent of k like in an s-wave superconductor and also there should be a normal-state-like 
component. These features are most naturally described by the three orbital band model ifTsf , where it is assumed 
that the superconductivity resides mainly in the 7 band, while the a and (3 bands may be considered in the normal 
O ■ state. Of course we cannot exclude the possible small superconducting order parameter associated with the a and /3 
band |i"9|j . In this perspective, the small magnetization seen by the muon spin resonance [^0| is rather puzzling. In 
addition, rather strong flux pinning in the vortex state is observed in Sr2Ru04 comparable to the one in £?-phase of 
UPt3 and Uo.97Tho.o3Bei3 pl[ . The latter two systems are considered to be in the nonunitary state p2| . Here it is 
worthwhile to mention that, although all these are consistent with the triplet pairing, both the magnetization and 
the strong vortex pinning should arise from some topological defects or dislocation in the samples. A recent small 
angle neutron scattering experiment in a magnetic field parallel to the c axis showed that the square vortex lattice 
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is almost everywhere in the B-T phase diagram |23|. This square vortex lattice is also very consistent with p-wave 
superconductivity we consider here p4]-p6|. 
In the following we take the order parameter 

A(k) = Ad(h ± ih) , (1) 

where d is assumed to be parallel to the c axis and k is the quasi-particle wave vector within the a-b plane. For 
simplicity, we also assume that the superconductivity resides only in the 7 band. This means the a and /3 bands 
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provide non-superconducting background. Recently starting from Eq.[l], microscopic studies of the vortex state p5| , p6| , 
the effect of impurities ^7j, and the quasi-particle spectrum around a single vortex p8[ have been discussed. 

Here we study spin dynamics in the p-wave superconductors with the order parameter given by EqJ3. It is important 
to realize that there are four distinct situations for the measurement of the dynamical spin susceptibility in two- 
dimensional p-wave superconductors. This is because the direction of the vectorial order parameter d(k) is fixed 
along a crystallographic axis which is perpendicular to the basal plane. One can apply a magnetic field along the 
directions perpendicular and parallel to the Ru02 plane. In each case, one can measure longitudinal and transverse 
susceptibilities. Therefore there exist four different susceptibilities. This is a unique property of the two dimensional 
p-wave superconductor because, in three dimensional case, the direction of the order parameter is not fixed and is 
always perpendicular to the applied magnetic field. 

We also assume that the pinning of d vector parallel to the c axis involves a finite pinning energy — ^XN&didz) 2 , 
most likely due to the spin-orbit coupling. Here we do not evaluate the exact value of £ld, but we expect that it will 
be rid(T) < A(T). Recently the pinning frequency was estimated by Tewordt. [^9| 

We believe that the experimental determination of Q d (T) will provide an insight in the pinning mechanism of d. In 
order to avoid future misunderstanding, we shall first explain the role of d here, d is called the spin vector which has 
been used in studying 3 He. It is perpendicular to the direction of the spin associated with the condensed pair 
Under these assumptions, we found the following results for the spin excitations. 

Case A: When magnetic field is parallel to the a-b plane (i.e., H _L d) 
We have the spin waves with the dispersion relations, 

«f = (1 - IfM + \f d {v F q) 2 (2) 

and 

u± = wjj + w| (3) 

for the longitudinal (u>\\) and the transverse resonance, respectively. The "pinning" frequency, £l d , is associated 
with the restoring tendency of the order parameter to the initial direction d against external perturbations. Here 
u>l = hbH is the Larmor frequency. I{— —jZo) is the dimensionless on-site Hubbard potential and f d is the dynamical 
superfluid density (lu ^ vfq) given by 

_ r ^A 2 tanh(#) 

fd -J A dE &7W=W> (4) 

where A is the magnitude of the full gap. These can be readily generalized in the presence of Fermi liquid corrections 
& 

Case B: When magnetic field is perpendicular to the a-b plane (i.e. H || d) 
The static susceptibility in this case is given by 

Y 

X = Xo 1 _ ly , (5) 

where Y = 1 — f s is the Yoshida function . 

There is no longitudinal collective mode, but a damped mode with the following dispersion exists in the longitudinal 
response. 

f 



UJ = i—v F \ q y{i-f d )i. (6) 

For the transverse response, we have the spin wave with the same dispersion relation as Eq|| 

The remainder of this paper is organized as follows. In section II, we show the calculation of the dynamical spin 
susceptibilities for four different cases mentioned above. We conclude in section III. Some of the technical details are 
relegated to the appendix. 
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II. DYNAMICAL SPIN SUSCEPTIBILITY AND SPIN WAVES 



Here we use Green's function method of Ref. M. The single particle Green's function in the Nambu space is given 



by 



G 1 (iuj n ,k) = iu n - 



P3 - A(fc • p)ai , 



(7) 



where pi and Ui are Pauli matrices acting on the particle-hole and spin spaces respectively, h is the unit vector parallel 
to the static magnetic field. u) n = (2n + l)nT is the fermionic Matsubara frequency, = ^l/(1 — I), £k = 5^7 — P, 
and A is the magnitude of the superconducting order parameter. Then the spin susceptibility is expressed as the auto- 



correlation function of spin operators. For example, the irreducible spin-spin correlation function, x 
can be computed from 



mi 



$(iu v , q) = I] Tr[a,G(p, uo n )a 3 G{ V - q, 



(8) 



where ot\ — P3<7i,ct2 — 02, <%3 — P3C3 are the spin vertices and come from the expression of spin density: Si ~ 
Trf^TQ/^], uj,, — 2vttT is the bosonic Matsubara frequency. 

In the fully renormalized spin susceptibilities, it is crucial to include the fluctuation of the superconducting order 
parameter related to the rotation of the d vector. Now let us look at the specific cases. 



A. Magnetic field is perpendicular to the order parameter: H || x 

1. Longitudinal susceptibility 

Since the direction of the magnetic field is in x-direction, the longitual susceptibility is given by Xxx ~ ([S X ,S X ]) 
and can be obtained from 

X° 

Xxx = 1 _ 7v° ' ^ 

A~xx 

where / is the exchange interaction and Xxx ^ s the susceptibility irreducible with respect to /. x xx includes the 
coupling between the spin and the order parameter, and can be written as 

= 00 , VxxgVxx (1Q) 

■\XX /^XX -1 TT ' \ J J 

1 - giixx 

where g is the strength of the interaction which is responsible for the superconductivity. Here Xx*L ~ {[Sx, S x ])o is 
the bare susceptibility, V xx ~ ([S Xl SA x ])q,V xx ~ ([5A X , S x ])o represent the coupling between the spin density and 
the order parameter fluctuation, and Tl xx ~ ([SA x ,SA x ])q the correlation between the order parameter fluctuations. 
In Nambu's notation, SA X ~ ATr[^T aip±a±^]. Notice that x xx consists of two parts - the quasiparticle contribution 
and the contribution from the excitation of the condensate. In the case of the longitudinal response, the £Il term in 
Eq.^ should be dropped out from the final expression. The details of the computation of each correlation function 
can be found in the appendix. We obtain (£ = vp • q) 

x £°>, q ) = iV(o)( c2 2 ~^ 2 / ) , 

V xx (u,q) = N(0)(^f) = V xat (u,q) , 

n ra ( W ,q)=.g- 1 -7V(0)(^^/) , (11) 

where 

f(u n = 4A 2 rr 2 - ^ f°° hf ^ {e/2T) {e ^ 4g2( " 2 ± e) ± 4C ' A2 (u) 

JK ,W ^ J a /P - A2 [(C 2 -w 2 ) 2 + 4£; 2 (a; 2 -C 2 )+4C 2 A 2 ] 2 -16w 2 £; 2 (C 2 -w 2 ) 2 ' V ' 
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and (F) = ^F(<f>). From now on, we set the density of state, N(0) = 1. In all these analysis we neglect 
Ql = wl/0- — I) for simplicity. 

Since the orientation of the order parameter is initially fixed along one of the crystallographic direction, there should 
be an energy scale, Qd, associated with the tendency to restore the original direction against external perturbations. 
This "pinning" frequency enters in the correlation between fluctuations of order parameters and leads to 

n xx {u,a) = g- 1 -{ e ~fj^ d f) ■ (13) 

Thus we get the following expression for X%xi 

X° = (V + ^ )) (14) 

In order to find a well-defined excitation, we consider the limit lu 3> v F q, where the quasiparticles do not generate 
dissipation which could make collective modes damped. In this case, the above expression can be simplified in the 
long wavelength limit (q <C k F ) as 



A XX 



i(v Fq ) 2 u 2 n 2 j d 



(15) 



2 LO 2 {\{v F q) 2 -uj 2 )(L0 2 ~\{v F q) 2 -^ 2 d ) 
Finally, in the limit u 3> v F q, we get the following full susceptibility after taking into account the exchange interaction 



u 2 £iy d + \{v F qf{u 2 -0 



2 ) 

d> (16) 



{\{v F q) 2 ~ u 2 )(lu 2 - \{v F qf - n 2 ) - I[±(v F q) 2 (iJ 2 - tl 2 ) + cu 2 n 2 f d } 
where fd = Um w ^oh m g^o(/) and h is given by 

fi= r dE ^m^. (17) 



\JE 2 - A 2 E 2 

It is called the "dynamical superfluid density" . Now we can read off the dispersion relation of a collective mode from 
the pole of the response function. We find that the longitudinal susceptibility supports a spin wave and its dispersion 
relation is given by 

UJ 2 = (l-If d )n 2 + ±fd(v F q) 2 . (18) 

This is consistent with the similar expression in superfluid 3 He-A phase @,|l|. Not surprisingly p-wave supercon- 
ductors have the longitudinal spin wave as in superfluid 3 He-A. 



2. Transverse susceptibility 

Now we analyze the transverse susceptibility. The transverse susceptibility, xa ~ ^(i^y + iS z , S y — iS z ]), can be 

obtained from 



(1 - ix°-+)x°+- + ix°++x°— 
(I - i x a +-)(i - ix\) - i 2 xi+x°- 



(19) 



where x?± ±\ are the susceptibilities irreducible to /. One can easily show that X+- = X-+ an( l X++ = X— — • Now 
X+- and X++ are given by 



A- = xT + E i 



V+kgVk- 



A7v 



k=y,z 

x° ++ = x° + \ + E > (20) 

, 1 — gUkk 

k=y,z 
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Herex+°_ ~ ([S v +iS z , S y -iS z ]) , X °++ ~ ([S y +iS z ,S y +iS z ]) are the bare susceptibilities. V ±k ~ ([S v ±iS z ,SA k }} , 
Vk± ~ ([SA k , S y ± iS z ])o and 11^^ ~ ([6Ak,6Ak])o are spin/order-parameter couplings and the fluctuation propagator 
of the order parameters respectively, where <5Afc ~ Trl^'akPicrt^]. 
When uol = 0, we find 



x°° + = 



C 



00 

A++ 



00 



V 



, 



II,. 



n, 



-n„ 



\ 4A 2 

In the presence of "pinning" frequency, Tl yy should be modified to 

iJ -w - .9 1 4A2 
Now incorporating finite and using Eq.EOL we get the following results 



(21) 



(22) 



Y° 



X- 



x° ++ = X°- = 



K Q 2 + n\ - u 2 ' 2 

(/) 1 ^ 2 (/) 2 

2 2 ((^_ C 2_ fi 2 )/) 



1 



^ 2 </> 5 



2<(^-C 2 -n 2 J/) 



(23) 



When lo 3> «f9, the transeverse susceptibilities irreducible to I in the long wavelength limit (q -C &f) become 



x+~ 



X- 



N(0) 
N(0) 



1 {vpg? 

2 uj 2 



u 2 m 2 



fd 



fd 



1 



{\{v F q)*-d*){d*-\(v F qy-V§ 2 2 w a - Q 2 - l(v F q) 2 



2 2Lu 2 -n 2 - \{v F q) 2 



(24) 



Using the above results and the RPA (random phase approximation) expression (see Eq.[19|) for the full transverse 
susceptibility, one finds two poles which correspond to a propagating spin wave mode and a damped mode. The 
dispersion relation of the spin wave is the same as the one in the longitudinal susceptibility: 



uj 2 = (l-If d )n 2 d + L0 2 +±f d (v F q) 2 , 



(25) 



where lol is the Larmor frequency. 



B. Magnetic field is parallel to the order parameter: H || z 

1. Longitudinal susceptibility 

If the magnetic field is applied along the direction of the order parameter, then the longitudinal susceptibility 
corresponds to \zz ~ ([S Z ,S Z ]}. As explained in previous sections, the full \ zz can be again obtained from 



A.ZZ 



i - ix° zz 



oo 

Azz 



l-gU z 



(26) 



z\)0, V zz 



Here x\ z 1S tne longitudinal susceptibility irreducible with respect to /. Notice that x zz ~ ([S Z ,S : 
{[S Z ,8A Z ]) ,V ZZ ~ {[6A Z ,S Z ]) , and U zz ~ ([SA Z , *A,]) 0j where SA Z ~ ATr[*t a3(0l(Tl ^]. 

In this geometry the superconducting order parameter does not move in the presence of the a.c. field. So there is 
no coupling between the spin density and the fluctuation of the order parameter. As a result, we have 



o 

A. z z 



1 - IX° Z z 



(27) 
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where 

,/Q _ ,,oo _ / C 2 (l ~ f) \ r< 

Xzz - Xzz - \—pi —) i Z ^) 

Eq. |27] gives a pole which gives the damped mode at 



LO = i-j=v F \qWl{\- f d ) . (29) 



2. Transverse susceptibility 

The transverse susceptibility, x~\ ~ \{[Sx + iS y ,S x — iS y ]), can be computed by using RPA expression (Eqjl9|) 

and the following relation 

00 , v +kgV k - 

— 



k—x,y 



gH-kk 



.o _ .,oo , y V+k9Vk± m 
±* i-gn kk K ' 



X++ = X++ 



Here X °°_ ~ {[S x + iS y ,S x - iS v }} , X + ° + ~ {[S„ + iS y ,S x + iS y ]) . Also V ±k ~ ([S x ± iS yi 6A k }} , V k± ~ 
([SAkjSx ± z5 !/ ])o and 11^^ ~ ([SA k ,SA k ])o are spin/order-parameter couplings and the fluctuation of the order 

parameter respectively. In Appendix B, we show that x+°-j V-{ i an d n_j are the same as Xxxi Vxx, an d rin when 

the magnetic field is along x. direction. It is also shown that x±± — V±± = IL±± = 0. Thus we can use the previous 
results, Eqjll] to get 



x Q + - = (^(c 2 + ., ;^ 02 )) (3D 



1 ( c 2 , 

y u; 2 -c 2 -n 2 

This is again exactly the same as Xxx> Eq.p^, when the magnetic field is along x direction. Since x+j 
the full susceptibility is just given by 



X°— 



X + - = • (32) 

This is also the same as the full longitudinal susceptibility Xxx for the case of H 1 1 x. p2| Therefore, x~\ m this case, 

supports a propagating spin wave mode (there is no damped mode) and the dispersion relation is given by in Eq, 18 
The transverse response X -\ when H || z is exactly the same as the longitudinal response for the case of H || x. 



III. CONCLUSION 

In this paper, we study dynamical spin susceptibilities in a two dimensional p-wave superconductor with A = 
d{k\ ± ik?) symmetry. This order parameter vector has been a strong candidate for the pairing symmetry of the 
superconducting Sr2Ru04. Due to the fact that the direction of the order parameter vector is fixed along a crys- 
tallographic direction, there are four possible susceptibility measurements: Longitudinal and transverse responses in 
the cases of the magnetic field parallel and perpendicular to the order parameter vector. The existence of spin wave 
modes in each case is examined and the dispersion relation is obtained. We found the following three modes. Two of 
them are spin waves. 

a) wf = (1 - If d )n} + y d (v F q) 2 , 

b) oj\ = LU 2 + uj\ , 

c ) u} = i^pklVW - fd) ■ 

The most crucial parameter is the pinning frequency. It is most likely that the spin wave is observable for electron 
spin resonance if -C A(T). We also believe that the experimental determination of f2^(T) will provide an important 
insight in the pinning mechanism of d. Naturally this will provide another test of p-wave superconductivity. 
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As it is pointed out elsewhere |3f], the superconductivity in Bechgard like (TMTSF) 2 PF 6 , (TMTSF) 2 C10 4 etc. 
is most likely of p-wave character as well. Therefore it is highly desirable to look for the spin wave in the above 
compounds as well. 

Unlike in superfluid 3 He, Sr2Ru04 is most likely in the vortex state. Otherwise the magnetic field would penetrate 
only in the surface. Nevertheless we believe that we can see the same expression for the spin wave dispersion, if we 
reinterpret the superfluid density with the one in the vortex state. We shall postpone the study of the spin wave in 
the vortex state to the future. 
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APPENDIX A: 

1. Magnetic field is perpendicular to the order parameter: H | x 

a. Longitudinal susceptibility 
Each correlation function can be computed from 

X° x °x(iuu,<l) =T^^Tr[aiG(p,w n )aiG(p-q,iw„-iuv)] , 

n p 

V xx (iuJv,q) = T^^Tr[aiG(p,w„)ai j oiCTiG(p - q, iw n - iuj v )\ , 

n p 

U xx (iuJ u , q) = T ^ X! Tr [aiPio-iG(p, w n )aipiC7iG(p - q, iu n - iw v )] . (Al) 

n p 

Summing over p (circular Fermi surface is assumed) leads to 



WtiWn+i/ + A 2 



xl° x (i^,q) = nTN(0)Y, ( 1 
V xx (iuj 1/ ,ci)=7TTN(0)Y^ 



^ + A 2 \M + , + A 2 ; (V^ + A 2 + v /c 2 +l/ + A 2 ) 2 + C 2 



V^ + A 2 Jc 2 +1 , + A 2 / V^+V^ + A 2 +C 2 



,,, x „ m x _ , , Un ^ v + A 2 \ <M + A 2 + y/<J> +v + A 2 

Iixx{tuv, q) = 7r7W(0) 2J 1 + ; , .2 1 ( A2 ) 

y/u*+A*Jul +v + A*J (v^+A^+J^ + A 2 ) +C 2 



where £ = ~vp ■ q and N(0) — mjlix is the two dimensional density of states. After summing over uj n and analytic 
continuation iu> v — > u> + iS, we get Eq.pT|. 

b. Transverse susceptibility 

Here x+°- ~ {[Sy + iS z ,S y — iS z ])o, X++ ~ {[S v + iS z ,S y + iS z ])o are the bare susceptibilities and they can be 
computed from 
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X°°_(iw„,q) = T^^Tr[a+G(p,w n )a-G(p-qJw n -iu; v )] , 

n p 

X^(iWv,q)=T^2lt[a + G(p,w n )a + G(p-q,iw n -iw I/ )] , (A3) 

n p 

where a± = -7^(0^2 ± io^s) respresenting the transverse spin vertices: -^={S V ± iS z ) ~ Tr^W*]. Similarly, we have 

V±fc(ia; w ,q) = T^^Tr[a ± G(p,w n )o; /i; piO'iG(p - q, iw n - iw„)] , 
n p 

Ukk(iu„,q) =T^2^2 r Tr[a k p 1 (7 1 G(p,uj n )a k p 1 a 1 G(p-q,iuj n -iuj 1/ )] . (A4) 

71 p 

It is also useful to represent these correlation functions in terms of Vt± t ±) ~ ([S y ± iS z ,SA±])o and II(±.±) ~ 
([6A±,6A±])o, where <5A± - Tr[^a±p 1 a l ^}. It is found that 

V±y = V y± = -j=(V± ± +V ±T ), 

V ±z = V z± = -^{±V ±± T V ±T ) , 

Kyy = ^(n++ + n+_ + n_+ + n__) , 

n« = -l(n ++ -n + _-n_ + + n__) , (A5) 

where 

^(±,±)(* w i',q) = ? 1 ^^Tr[a±G(p,a;„)a± / 9iCTiG(p - q, iw„ - iav)] , 

n p 

II(±,±)(^,q) = T } } r Tr[a±piaiG(p,ui n )a±piaiG(p - g,,ioJ n - iwu)] ■ (A6) 

n p 

After summing over p, we get 



v2 



xr + =x--=*™(o)£ 



-A 2 \ V^n + A2 + V ^n+„ + A 



V^ + A 2 ^n + , + AV (v^TA^+^^ + A 2 ) +C 2 
V+_ = = + = = TrTiV(O) ^ 



1 -iuj u A \ yfi% + A 2 + + A 

2 



n+- = n_+ = ^nv(o)^ 



A 2 \ V^n + A 2 + J <4 +!/ + A 



V / ^+A 2 ^ 2 +,+A 2 ; ( V / W 2 +A 2 + / u ,2 +y + A 2') +C 2 



2 



n ++ = E__ = TrTiV(O) £ 1 + " ^ = = v , (A7) 

n V V^+A^V W "+- + A / (v / ^ 7 A 2 + V^+- + A 2 ) +C 
Now summing over and analytic continuation iuj nu — > w + iS lead to Eq.|2l|. 

APPENDIX B: 
1. Magnetic field is parallel to the order parameter: H || z 
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a. Longitudinal susceptibility 



Notice that X ° z ° z ~ {[S Z ,S Z ]) , V zz ~ {[S Z ,SA Z ]) ,V ZZ ~ {[6A Z ,S Z ]) , and Tl zz ~ ([5A 2 , 0A,]) O , where <5A Z ~ 
ATr[^E , ^o;3pi(Ti v E']. These correlations functions can be expressed as 

X°°(*w V) q) =T^^Tr[a 3 G(p,w„)a 3 G(p-q,iw„-iw I ,)] , 

n p 

= T^^Tr[a 3 G(p 

n p 

n zz (Mj„,q) = T^^Tr[a 3 picriG(p,w n )Q;3picriG(p - q, iu n - . (Bl) 



^From the consideration of the matrix elements, we can easily see that V zz = V zz = 0. Therefore, there is no mixing 
between spin fluctuations and the order paramter in this case and x° zz is just given by x zz . After summing over p, 
we get 



oo 

A. z z 



{iu u , q) = irTN(0) 2_J 1 7= , . .2 • ( B2 ) 



Summing over oj„ and analytic continuation lead to Eq, 28 



fe. Transverse susceptibility 

Here ~ ([S x + iS y ,S x - iS y ]) Q , X ++ ~ ([S x + iS y ,S x + iS y }) . Also V ±k ~ ([5 X ± iS B , <JA fc ]) , T4± ~ 
[SAk, S x ± i5 B ])o, and likk ~ ([<5A/j, <5Afc])o. These correlation function can be again calculated from 

X %(iuv, <1)=TJ2Y1 T ^+ G ^ "n)P T G(p - q, icu n - iu v )} , 

n p 

V±fc(i^,q) = r^^Tr[/3±G(p,w„)o!fc/9icriG(p- q, iw„ - iu v )] , 

n p 

Ilfcfc^cJv^) = Ty^y^Tr[akpi(TiG(p,u) n )akPi(TiG(p - q,iu> n - iu> v )] , (B3) 



where /3± — -^{ a i ± ^2) respresenting the transverse spin vertices: -^{S x ± iS^) ~ Trf^T^^r]. One can also 
rewrite these correlation functions in terms of V^± ±^ ~ {[S x ± iS^, 5A±])o n( ±i± ) ~ ([SA±, SA±])q, where SA± ~ 

^ = F„± = -^(v±± + V± T ) , 

V±y = V z± = -L(±V ±±T V ±T ) , 

V 2« 

n ax = i(n ++ + n+_ + n_+ + n__) , 

n ra = -i(n ++ -n+_ -n_ + + n__) , (B4) 

V(±,±)(iu) u ,q) =T^^Tr[^ ± G(p,cj„)/3±p 1 (T 1 G(p-q,iw n -icj Iy )] , 
n p 

n (±,±)(«^,q) = T^^Tr[/3 ± pi(TiG(p, w„)/3±picriG(p - q,iw„ - iu v )] . (B5) 



where 



n p 

From the matrix elements, one can see that 

,00 _ ,00 



x7+ = x { -- = o, v ++ = V— = o, n ++ = n__ = o. (B6) 
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Thus Eq.BfJ can be simplified as 



X°- 



X° + °- 



V+-9V+- 

i - «?n + _ 



(B7) 



After summing over p, we obtain the same equation as Eq. A2. Thus we can use the previous results (EqJTl]) to get 
Eq.[n] which is same as Eq@ 
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